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In superconductors, the search for special vortex states such as giant vortices focuses on laterally 
confined or nanopatterned thin superconducting films, disks, rings, or polygons. We examine the 
possibility to realize giant vortex states and states with non-uniform vorticity on a superconducting 
spherical nanoshell, due to the interplay of the topology and the applied magnetic field. We derive 
the phase diagram and identify where, as a function of the applied magnetic field, the shell thickness 
and the shell radius, these different vortex phases occur. Moreover, the curved geometry allows 



these states (or a vortex lattice) to coexist with a Meissner state, on the same curved film. We 
2 have examined the dynamics of the decay of giant vortices or states with non-uniform vorticity 

into a vortex lattice, when the magnetic field is adapted so that a phase boundary is crossed. 
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I. INTRODUCTION 



Quantized vortices are a quintessential property of superfluids and superconductors. The 
energetically favored state when multiple quanta of vorticity are present, is a lattice of 
singly quantized vortices. In ultracold Fermi gases, the recent observation of such a vortex 
lattice formed the 'smoking gun' proof for superfluidity [Ij. In nanoscopic superconducting 
samples, controlling the vortex behavior is essential for the development of new devices based 
on fiuxon dynamics [2] . The confinement of Cooper pairs to the length scales comparable to 
the correlation length also offers the prospect to probe fundamentally new phase topologies 
predicted by the theory, such as giant 0, H] and ring-like vortices [S] . This has led to renewed 
experimental efforts to observe giant vortex states, both in superconductors [HI [7] and in 
superfiuid atomic gases |8]. 

In this contribution, we argue that superconducting spherical nanoshells form a promis- 
ing candidate for realizing giant vortex states, and for engineering phase transitions between 
those states and a vortex lattice. Moreover, we show that nanoshells allow the co-existence 
of a Meissner state and a vortex state in equilibrium on one and the same superconducting 
film. Nanoshells are hybrid nanostructures consisting of a dielectric core (usually a silicon 
oxide nanograin), coated with a thin layer of metal[9]. When the metal in its bulk form is a 
superconductor, the nanoshell below the critical temperature will also exhibit superconduc- 
tivity in the thin shell around the isolating core. 

The superconducting order parameter in the nanoshell is well described by a macroscopic 
wave function tp = \ip \ e**^ that obeys the coupled time-dependent Ginzburg-Landau (TDGL) 
equations. Vortices are characterized as topological defects in the phase (requiring a 
vanishing gap For thin shells, the description is simplified in two important ways. 

Firstly, when the shell thickness is much smaller than the London penetration depth, the 
magnetic field will be only weakly perturbed by the nanoshell. Secondly, when the shell is 
thinner than the coherence length, the order parameter ip will not vary substantially in the 
radial direction in the shell; that is, ip will only depend on the spherical angles Vt = {6, 0}. 
In the radial direction, ip will be constant in the shell, and zero outside it. Note that 
confining ip to the shell leads to an effective Ginzburg-Landau parameter n that differs from 
its bulk value. In Section II we present the formalism, and in Sec. Ill the results, for thin 
shells. When the shell thickness is increased and becomes non-negligible with respect to 
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the penetration depth, the magnetic field will be more strongly perturbed, and the field 
gradients affect the energetics. This case and the effect on the phase diagram are discussed 
in Sec. IV. Finally, we summarize the results for vortices in nanoshells in Sec. V. 

II. GINZBURG-LANDAU FORMALISM ON THIN SHELLS 

We assume that the shell is sufficiently thin for neglecting variations of the order parame- 
ter across the shell. In other words, the order parameter ip will only depend on the spherical 
angles = {6,(j)}. We use the spherical coordinates r, 6, with the origin at the center 
of the sphere. The angle 6 is counted from the z-axis parallel to the external homogeneous 
magnetic field. Like in Ref. [121 ^6 will make the used variables dimensionless by expressing 
lengths in units of V2X, magnetic fields in units of ^q/ (AnX^), and the vector potential in 
units of $o/(2"\/2vrA), where A is the penetration depth, $o = h/{2e) is the magnetic fiux 
quantum, h is the Planck constant, and e is the elementary charge. Thus, the dimensionless 
parameters W, and H are linked to the radius of the nanoshell TZ, its thickness W, 
and the applied magnetic field H by the expressions R = Tl/{\/2X), W = >V/(-\/2A), and 
H = 47rA^7i/$o, respectively. 

In our numerical treatment of superconducting states on spherical shells we exploit the 
time- dependent Ginzburg-Landau equation, which is known to be a powerful tool for study- 
ing both the dynamic and static properties of superconductors. For a thin shell under 
consideration, the behavior of the order parameter in a fixed (or slowly varying) magnetic 
field can be described by the TDGL equation (cp. [TOl ITT] ) 



where k is the Ginzburg-Landau parameter, A is the (dimensionless) vector potential, and 
= eg{d/d9) + sm^^{6){d/d(j)). The dimensionless variable r is linked to the time t by 
the relation r = Dt/TZ"^, with D, the normal-state diffusion constant. 

The vector potential A can be represented as a sum of the contribution Ai, related 
to supercurrents in the shell, and the contribution Aq, which corresponds to the external 
magnetic field H. The vector potential Aq is chosen in the form 



dijj 



(Vn - iRP^f ^ + 2(fi:i?) V(l - IV^n 



(1) 



Ao = 



HrsinO 



(2) 



2 



3 



In the case of a constant applied magnetic field H, with increasing r the function V', given 
by Eq. ([T]), approaches one of the (meta)stable states of the system {dip/dr 0). The 
thermodynamically stable state is to be found by comparing the Gibbs free energy for 
different solutions. The difference in the Gibbs free energy between a superconducting state 
and the normal state at the same magnetic field is given by the equation 

27V n 

^G = ^^ j d<p j de[A,.]- /c^iv^r] , (3) 



where AGq corresponds to the superconducting state with no vortices at if = 0, i.e. the 
Meissner state present on the complete surface. The dimensionless density of supercurrents 
is denoted by j and expressed in units of <l>oc/(8-\/27r^A'^). 

In this section, the shell is assumed to be sufficiently thin in order to make negligible the 
magnetic fields, induced by supercurrents. Correspondingly, we can neglect Ai as compared 
to Aq. Then, as seen from Eqs. ([T| and (|2|, two independent parameters, which govern the 
solution of Eq. ([T]), remain: 

• the dimensionless size of the nanoshell p = kR = TZ/(\^^), determined by the ratio 
of the shell radius TZ to the Ginzburg-Landau coherence length ^, and 

• the parameter rj = = -kI-CR? /^q, equal to the number of flux quanta of the 
applied field that pass through the equatorial plane of the sphere. 

When the magnetic field is increased beyond a critical value (computed below), a first 
vortex appears for nanospheres with radius large enough to sustain the vortex core, as 
depicted in Fig. [l| panel (a). Upon further increasing the magnetic field, more quanta of 
flux can penetrates the spherical surface. This can be accomodated in a variety of ways: 
for example as a giant vortex carrying more than one quantum <l>0; shown in panel (b) of 
Fig. [TJ In this case, the angular momentum is uniform over the spherical surface. It is also 
possible to envisage states with non-uniform distributions of angular momentum: a value 
Li near the poles and a value L2 in a band near the equator. Such states are characterized 
by a ring-like vortex separating the regions with different angular momentum, as illustrated 
in panel (c) of Fig. [T| Also states that do not have axial symmetry should be investigated: 
we will show that these are in many case the stablest state and that they then consist of 
an array of singly quantized vortices, illustrated in panel (d) of Fig. [l] Such states will be 
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single cD„-vortex ring-like vortex 




FIG. 1: A few examples of vortex structures on a spherical shell are illustrated in this figure, to 
clarify the nomenclature used in the text. The amplitude and phase of the order parameter are 
shown as the saturation and hue of the color scale, respectively, in such a way that the vortex 
core region is white. Panel (a) shows a single <&o-vortex. In panel (b) a giant vortex (carrying 
multiple quanta of flux) is depicted. Panel (c) illustates a ring-like vortex separating regions with 
different angular momentum. Finally, panel (d) shows the <I>o-™ultivortex state, where an array of 
singly-quantized vortices is present. 

denoted by $o-iiiultivortex states, to emphasize that every vortex carries a single quantum 
of flux. In the next subsections we start by investigating the axially symmetric states: giant 
vortices and ring-like vortices. Then, in the next section we investigate the condition under 
which those states decay into $o-niultivortex states, and the dynamics of this decay. 



A. Giant vortex states 

First, let us consider superconducting states, which keep the axial symmetry of the sys- 
tem, so that the order parameter ip can be written in the form ip = /(^)exp(iL0). where L 
has the sense of the winding number (vorticity). Then for a stationary distribution f{6) the 
Ginzburg-Landau equation ([T| reduces to the one-dimensional equation 
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FIG. 2: The free energy difference is shown as a function of r/, the number of flux quanta that pass 
through the equatorial plane of the sphere with radius p = 8, for axially symmetric states. The 
black solid curves are for states with uniform vorticity (and a giant vortex). The red dotted lines 
are for states with ring-like vortices. 

with boundary conditions, determined by the requirement that the ^-component of the cur- 
rent density must be zero at the 2;-axis: df /dO]^^^^ = 0. SoUd lines in Fig. |2] illustrate 
typical behavior of the free-energy difference AG as a function of t] for cylindrically sym- 
metric states with different vorticity L. In the case of a thin spherical shell with p = 8, 
as illustrated in panel (a) of Fig. [2| the value of L in the lowest cylindrically symmetric 
state increases with rj from at = to 9 at = 12.5. The modulus / and phase L(j) are 
illustated in panel (b) of Fig. [l| using hue and saturation of the color scale respectively. An 
increase of the applied magnetic field is seen to result also in a significant increase of the 
Gibbs free energy of the lowest state. 

B. Ring-like vortices 

In Ref. [12], when analyzing supercondicting states in hollow cylinders, it was suggested 
that - under certain conditions - cylindrically symmetrical states with changing winding 
number can be more energetically favorable than the states with uniform L. Our calculations 
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show that a similar situation occurs also in thin spherical shells with the dimensionless size 
larger than p ^ 6 (i.e., for TZ > 8.5^), but as we will show in the next section, such states 
decay into a lattice of singly-quantized vortices breaking the cylindrical symmetry. 

We have compared the Gibbs free energies for axially symmetric states with uniform 
winding number L and those for states where Li, the winding number at < ^ < and 
7T — 00 < 9 < vr, differs from L2, the winding number at 60 < 6 < n — Oq. The order 
parameter of the latter states on the sphere are illustrated in panel (c) of Fig. 1 - we will 
refer to such states as "ring-like vortex" states. The continuity of the order parameter as 
a function of 6 requires vanishing f{R,6) at the boundaries between regions with different 
winding numbers, i.e., at 6 = 9q and 6 = n — 60, as can be seen in panel (c) of Fig. [!} 
At sufficiently strong magnetic fields the Gibbs free energy for states with ring-like vortices 
can become lower than that for states characterized by a unique winding number L over 
the whole 6 range. This is illustrated by Fig. [2| where the dotted lines show the calculated 
free-energy difference AG{ri) for states with Li = 1 and different L2 on a shell with p = 8 in 
the case of 6q = O.ITStt. It is worth mentioning that the values ^q™™'', which minimize AG 
for the lowest ring-like vortex state at a given 77, are rather insensitive to the dimensionless 
nanoshell size p (at least, for p < 10). At the same time, the parameter is an increasing 
function of rj. Thus, our calculations show that this parameter changes from ^q™"'' ^ 0.127r 
to ^ 0.27r when increasing rj from 7 to 15. However, moderate variations of Oq around 

^Q™™^ only slightly affect AG for the lowest state with ring-like vortices. That is why in 
Fig. |2]we restricted ourselves to the case of a fixed value Oq = 0.1757r, which coincides with 
at r/ = 10. 

In Fig. |2j the curve labeled with '1,1' corresponds to the state with a ring-like vortex and 
uniform vorticity, which is qualitativley similar to the states analyzed in Ref. |5l The free 
energy of this state is always significantly higher than AG for the lowest giant vortex states 
with point-like core. As can be further seen from Fig. [2| at r/ > 10 the ring-like vortex states 
(Li 7^ L2) appear the most energetically favorable among the axially symmetric states. 
Thus, at ?7 > 12.5, the difference in AG/AGq between the state with L = 9 and the state 
with Li = 1, L2 = 9 is larger than 0.025. At even larger values of rj, states with three 
different regions of vorticity Li, 1/2,-^3 (characterized by two ring-like vortices) can become 
stable. However, our calculations show that in shells with 1] > Q, where such ring-like 
vortices allow for decreasing the free energy of giant vortex states as compared to the case 
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of a giant vortex, even lower values of AG can be achieved by breaking up the ring like 
vortex (or vortices) into an array of singly-quantized vortices. A natural question arises of 
how stable are the aforedescribed giant vortex states with respect to decay into multiple 
singly quantized vortices. In order to answer this question, one has to return to the TDGL 
equation ([T]). 

C. Numerical treatment 

The finite-difference scheme, applied here to solve Eq. ([T]), is similar to that of Ref. [H], 
with necessary adaptations to the case of a spherical 2D-system. Two-dimensional grids, 
used in our calculations, typically have > 100 equally spaced nodes in the ^-interval from 
to 7r and > 150 equally spaced nodes in the 0-interval from to 27i. Cyclic boundary 
conditions for ip are applied at = and = 27r. The boundary conditions at ^ = 
and 6 = n are determined by the requirement ip\g^Q^ = const (0). The step of the time 
variable r is automatically adapted in the course of calculation. This adaptation is aimed to 
minimize the number of steps in r, necessary for approaching a steady solution of Eq. ([T]), 
and - at the same time - to keep the solving procedure convergent. On average, the step in 
r is ~ 10^^ to ~ 10""^ depending on the used grid as well as on p and t]. When starting at 
r = from a random distribution of ip (with \-ip\ <^ 1), a (meta)stable solution of Eq. ([T]) 
is achieved typically at r < 100. When analyzing (meta) stability of states in a spherical 
shell, one has to keep in mind that a transition between states with different vorticity, in 
general, requires symmetry breaking. This means that simulations, which assume a perfectly 
symmetric spherical nanoshell, would tend to overestimate stability of a state with respect 
to a possible transition to another state with lower free energy. In order to model the effect 
of imperfections, inevitably present in realistic nanoshells, we consider spherical shells with 
small angular variations 6p{6, (p) of the parameter p. Importantly, for relative magnitudes 
\Sp\/p ranging roughly from ~ 10^*^ to ~ 10^'^, the results of simulations practically do not 
depend on a specific choice of the magnitude and distribution of these non-homogeneities. 
An appreciable effect of those imperfections on stable distributions of the order parameter 
appears only for \Sp\/p > 0.1. 
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FIG. 3: Evolution of the angular distribution of the squared modulus of the order parameter in 
a thin spherical superconducting shell with p = 8 at = 10 in the case when the initial state (at 
r = 0) is a giant vortex with L = 7. Different panels correspond to different time r. 

III. RESULTS AND DISCUSSION FOR THIN SHELLS 
A. Decay of giant and ring-like vortices 



In order to examine the stability of giant vortex states with respect to decay into multi- 
vortex states, we apply the computation scheme described in the previous section, starting 
at r = from a distribution of that corresponds to a giant or ring-like vortex state. 
Typical examples of the evolution of the order parameter distributions are shown in Figs. [3] 
and |4] for the cases when the initial state is a giant vortex and a ring-like vortex, respec- 
tively. In the case of p = 8 and 1] = 10, the thermodynamically stable state corresponds to 
7 pairs of vortices with a single quantum $o of flux each, and it has a relative free energy 
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AG/AGo ~ —0.812, approximately 0.07 lower than the value of AG/AGq for the lowest 
giant vortex state (see Fig. |2]). Such lattices of vortices with each a singly flux quantum $0 
will be denoted as "<l>o-niultivortex states" . 

As illustrated by Figs. [3|a) to|3](c), within a r- interval ~ 1 the initial giant vortex state 
with L = 7 transforms into a chain of 7 singly quantized vortices, which surround each pole 
of the sphere. In the course of the further rearrangement of the vortex pattern, one of the 
vortices moves to the pole, while the remaining 6 vortices tend to form a symmetric chain 
around the pole (see Figs, [sj^d) to[3]^f)). A free-energy gain due to this rearrangement is by 
one order of magnitude smaller than that due to the decay of the initial giant vortex into 
single vortices. Correspondingly, the r-interval, necessary for this rearrangement, appears to 
be relatively long: only at r ^ 10 the solution reaches the equilibrium symmetric configura- 
tion of vortices (not shown in Fig. |3]), similar to that found in Ref. |13l where $o-niultivortex 
states on a thin hollow sphere were studied in detail. As seen from Fig. |4| the transition 
from a ring-like vortex state (Li = 1 for < 9/n < 0.175 and 0.825 < 9/n < 1; = 7 
for 0.175 < 6/tt < 0.825) to a $o-niultivortex state is even faster. The ring-like vortex core, 
which is present in the initial state [see Fig. |4]^a)], decays into a chain of 6 single vortices 
very quickly: clear signatures of this decay can be found already at r < 0.01 [see Fig. |4Fb)]. 




The equilibrium state with a vortex at the pole and 6 vortices, symmetrically surrounding 
the pole, is formed already at r = 0.2. 

The results of our calculations clearly indicate that giant and ring-like vortex states are 
rather unstable in spherical shells with relatively large p. This does not mean, however, 
that giant vortex states on a spherical shell are never stable. A decrease of the shell radius 
and/or an increase of the applied magnetic field enhance the role of the Lorenz forces, which 
act on the supercurrents and tend to drive vortices towards the poles of the shell. As a 
result, for sufficiently small p and sufficiently large 77, the distance between vortex cores 
in a $o-iiiultivortex state becomes so small, that physically a $o-iiiultivortex state appears 
undistinguishable from the corresponding giant vortex state. A similar continuous transition 
from a $o-niultivortex state to a giant vortex states with increasing magnetic field was 
recently found when solving the linearized Ginzburg-Landau equation for superconducting 
spherical grains [14j. Of course, in the case of such a continuous transition, the boundary 
between thermodynamically stable $o-iiiultivortex states and giant vortex states can be 
drawn only approximately. As a criterion of a transition from a multivortex state to a giant 
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FIG. 4: Evolution of the angular distribution of the squared modulus of the order parameter in 
a thin spherical superconducting shell with p = 8 at = 10 in the case when the initial state (at 
r = 0) is a ring-like vortex state with Li = 1 for < j-K < 0.175 and 0.825 < 9/tt < 1; L2 = 7 for 
0.175 < O/tt < 0.825. Different panels correspond to different time r. 



vortex state, here we have chosen the condition that the angular distance of vortex cores 
from the pole becomes smaller than (lOp)^^. 



B. Phase diagram for thin shells 

Our results, related to thermodynamically stable states on thin spherical shells, are sum- 
marized in Fig. [5| where the solid lines indicate boundaries of stability regions for the 
normal state, the superconducting Meissner states, single $o-vortex states, giant vortex and 
•I'o-niultivortex states. The dashed line indicates the boundary between the regions, where 
giant vortex states (to the left from this line) or $o-Hiultivortex states (to the right from 
this line) are thermodynamically stable. As seen from Fig. [sj formation of vortices can be 
energetically advantageous only on sufficiently large shells: at p > 0.63 for states with L = 1, 
at p > 0.85 for giant vortex states, and at p > 1.95 for $o-niultivortex states. 

While in Fig. |5] the phase diagram for thin spherical superconducting shells is shown in 
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FIG. 5: Phase diagram for thin spherical superconducting shells in the (p, ry)-plane. The boundaries 
between the regions with the thermodynamically stable normal state, the Meissner state, the single 
$o-vortex state, giant vortex state, and <I>o-™ultivortex states are shown by solid lines. The dashed 
line approximately indicates the boundary between the regions, where giant vortex states or ^q- 
multivortex states are the thermodynamically stable state. 

the (p, r])-plane, it seems interesting to analyze the phase boundaries also in more common 
form: in terms of the applied magnetic field Ti. and the temperature T. We assume that 
the temperature dependence of the penetration depth A is described by the empirical rela- 
tion A(T) = A(0)/ a/1 — (T/Tc)^, while the (less important) temperature dependence of the 
Ginzburg-Landau parameter k, is roughly given by the expression k(T) = k(0)/[1 + (T/Tc)^] 
(see, e.g., Ref. [TSil . In Fig. |6] we plot the phase boundaries for thermodynamically stable 
normal states, Meissner states, single $o-vortex, giant vortex, and $o-multivortex states 
on spherical superconducting shells with different radius TZ, measured in units of the zero- 
temperature Ginzburg-Landau coherence length ^(0). As illustrated in Fig. [6|^a), in the 
case, where the radius TZ is much larger than ^(0), giant vortex states are thermodynam- 
ically stable (at moderate applied magnetic fields) only in the close vicinity of the critical 
temperature Tc. With decreasing TZ, the stability range of giant vortex states gradually 
extends towards lower temperatures and lower values of the magnetic flux through the shell. 
Correspondingly, in sufficiently small shells the stability range of $o-multivortex states is 
restricted to a relatively narrow interval of Ti and to T, significantly lower than [see 
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FIG. 6: Phase boundaries for thin spherical superconducting layers with different radius 7^ as a 
function of the temperature T and the applied magnetic field H. Boundaries between the thermo- 
dynamically stable normal state, Meissner state, single $o-vortex state, giant vortex state, and $o- 
multivortex states are shown by solid lines. The dashed line approximately indicates the boundary 
between the regions, where giant vortex states or ^o-multivortex states are the thermodynamically 
stable state. 
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Fig. [6[b)]. On even smaller shells (with TZ close to ^(0)) no stable $o-niultivortex states are 
possible [see Fig. |6](c)]. For those small shells, the superconducting phase persists only for 
relatively weak (few $o) magnetic fluxes through the shell. At the same time, as implied by 
a comparison between panels (b) and (c) in Fig. [6| the values of the applied magnetic field 
7i, which correspond to transitions with an increase of vorticity by 1, become significantly 
higher when decreasing the shell size. 



IV. VORTEX STATES ON THICK SHELLS 
A. Magnetization effects in thick shells 

Now, let us extend our analysis to the case of relatively thick spherical shells, where 
magnetic fields, induced by supercurrents, which flow in a shell, are of non-negligible. At 
the same time, we assume that the thickness of a shell is still sufficiently small for neglecting 
variations of the order parameter and of the vector potential A = Aq + Ai across the 
layer. For such a shell, also currents across the layer can be neglected. Expressing the vector 
potential Ai through the density of current j as 

Ai(r) = — / rfV^^. (5) 
27r y |r - r'l ^ ' 

the non-negligible components of the product i?Ai, which enters Eq. ([T]), can be written 
down in the following form: 

7r 2-77 

RAie = / de' sin 9' [ dcj)' { [sin ^sin^' + cos 9 cos 9' cos (0 - 0')] Rje' {9', (j)') 

2v27r J J 



-|- cos 9 sin (0 — (f)') Rj^r {9' , </)')} [1 — cos 9 cos 9' — sin 9 sin9' cos (0 — 0')] ,(6) 



TT 2tT 



RAi^ = — ^— / d9' sin 9' / d(f)' {cos 9' sin {(p' — (f)) Rje 
2v27r J J 



+ cos (0 - 0') Rj^' {9', 0')} [1 - cos 9 cos 9' - sin 9 sin 9' cos (0 - 0')]"^^^ , (7) 



where W is the dimensionless thickness of the shell. On the other hand, in the case of 
constant or slowly varying magnetic fields, using the relation 



j = Re 



(8) 
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the products Rjg and -Rj^, which enter Eqs. ^ and ([7|, can be expressed through ijj, RAg, 
and RA^i, as 



Rjg = Im 



^ de 



RA^g\ij\\ (9) 



sm & 



-{t] sine + RA,^)\^Ij\^. (10) 



In order to find the order parameter and the corresponding vector potential, we solve 
self-consistently the set of equations ([T]), ([6]), and ([T]), using relations (|9|, and (10). From 
Eqs. ([T]), (|6]), ([t]), (|9]), and ( 10 ), one can see that for relatively thick shells under consideration 
a set of independent parameters, which govern the solution, can be chosen as r], p, and u, 
where the introduced additional parameter u = WR = W7?./(2A^) is linearly proportional 
to the thickness of the nanoshell and to its radius. 

Figure [7] gives few examples of magnetic- field distributions, which correspond to thermo- 
dynamically stable states in spherical shells with p = 8, a; = 30 and p = 8, uj = 10. The 
patterns of magnetic-field lines, displayed in Fig. [7[ are plotted for the particular case of the 
Ginzburg-Landau parameter k = 0.8, the (mean) dimensionless radius of the shell R = 10, 
and the dimensionless thickness W = 3 and W = 1. In general, none of the three mutually 
orthogonal components of the magnetic field B = V x A is zero, so that the field lines are 
"three-dimensional" . In Fig. [7| however, we restrict ourselves to field-line patterns within 
symmetry planes, where the field lines are "fiat". As seen from Fig. [7]^a), even in the case 
of a relatively thick shell {W = 3) the magnetic fields, induced by the supercurrents in the 
Meissner state, are not sufficient for complete screening of the applied magnetic field inside 
the shell. Nevertheless, not only in the case of W = 3 but also for a significantly thinner 
shell with W = 1 [Fig. [7|^b)], the net field inside the shell is much weaker than H. In the 
case of the state with L = 1, the magnetic fiux, captured by a vortex pair, is seen as an 
increased density of field lines at the poles of the sphere [Fig. [7]^c)]. At the same time, in 
the depth of the sphere the magnetic field is relatively homogeneous, only slightly increasing 
towards the 2;-axis. Also for states with higher vorticity, a considerable local increase of the 
magnetic-fiux density takes place only at the vortex cores within the superconducting shell, 
while in the depth of the sphere the density of magnetic-field lines is considerably more 
homogeneous [see panels (d) to (f) in Fig. [?]. 
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FIG. 7: Magnetic-field lines and distribution of the squared modulus of the order parameter for 
thermodynamically stable states in superconducting spherical shells with k = 0.8, R = 10, W = 3 
[panel (a)] and k = 0.8, R = 10, W = 1 [panels (b) to (f)]. The results are shown for the xz-cross- 
section [panels (a), (c), (d), (f)] and the yz-cross-section [panels (b), (e)] of the shell for different 
values of the parameter r] and vorticity L: rj = 5, L = [panels (a) and (b)], r] = 6.5, L = 1 [panel 
(c)], r] = 8, L = 2 [panels (d) and (e)], rj = 9.5, L = 3 [panel (f)]. Insets: angular distributions 
of the squared modulus of the order parameter for the same values of t] and other parameters. 
Vertical dashed lines on each inset correspond to the cross-section displayed on the main panel. 

B. Phase diagram for thick shells 



As seen from Fig. [7| the magnetic fields, induced by supercurrents, can be considerably 
large even for shells with quite moderate thickness {W ~ 1). These fields strongly affect 
the stability range for superconducting states with different vorticity in a spherical shell. 
In Fig. |8j we present the calculated phase diagram for relatively thick spherical shells with 
o; = 10. As follows from a comparison of Fig. |8] to Fig. [5} an increase of the thickness 
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FIG. 8: Phase diagram for thick (u = 10) spherical superconducting shells in the (p, r/)-plane. 
The boundaries between the thermodynamically stable normal state, Meissner state, single $o- 
vortex state, giant vortex state, and <I>o-™ultivortex states are shown by solid lines. The dashed 
line approximately indicates the boundary between the regions, where giant vortex states or ^q- 
multivortex states are the thermodynamically stable state. 

of a spherical shell results in a well-pronounced shift of the boundaries between states with 
different vorticity towards higher magnetic fields rj. In particular, for p > 0.8, the range of 
rj, where Meissner states are thermodynamically stable, is more than two times wider in the 
case of = 10 as compared to the case of u 0. One can also see that for a relatively 
thick spherical shell {uj = 10) the boundary between giant vortex and $o-Kiultivortex states 
is shifted towards significantly larger values of p as compared to those in the case of — 0. 
The increased stability of giant vortex states agree with the results recently obtained by 
Baelus et al. [H] for the limit W ^ R oi a. full sphere, in the framework of linearized 
Ginzburg-Landau equations. 

In Fig. [9] the phase boundaries for thermodynamically stable normal states, Meissner 
states, single $o-vortex, giant vortex, and $o-niultivortex state states are plotted in the 
(T, 7i)-plane. The results are shown for shells with different radius TZ and thickness W. In 
order to keep the plots more universal, it is convenient to express the thickness in units of 
.^(0)/€^(0). For thick nanoshells, there is a much more pronounced increase of the transition 
fields, which correspond to a change of vorticity, with lowering temperature [cp. Figs. [9|^a) 
and |9](c) to Figs. |6]^a) and[9]^c)]. When comparing Fig. [oj^a) to Fig. [6]^a) one can also see 
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that with increasing the nanoshell thickness the temperature range, where giant vortices are 
thermodynamically stable, extends towards lower temperatures. With decreasing the shell 
radius, this effect becomes quite pronounced even for relatively small values of W/[^(0)k^(0)] 
[cp. Fig. [9|b) to Fig. [6|b)]. In sufficiently thick nanoshells, the temperature range, where 
$o-niultivortex states are thermodynamically stable, reduces to zero [see Fig. [9|^c)], although 
in thin nanoshells of the same radius this range is relatively wide [see Fig. [6]^b)]. Of course, 
when the temperature approaches T^, the phase boundaries become almost insensitive to 
the value of W. Indeed, at T Tc the parameter uj always goes to zero [due to an increase 
of the penetration depth A(T)], so that any nanoshell appears effectively thin. 



V. CONCLUSIONS 



Curving a superconducting film into a spherical shell changes its vortex-related properties 
drastically due to topological constraints. The hairy-sphere theorem |T6] is a straightforward 
example of such a constraint: it states that, in contrast to the situation on a flat film, there 
exists no nonvanishing continuous tangent vector field on the sphere. So, every nonvanishing 
supercurrent velocity field requires discontinuities, such as vortices. The interplay between 
the Lorentz force due to an applied field and the vortex superflow will force these vortices 
away from the equator (leaving an equatorial "Meissner band") and towards the poles. This 
results in a 'polar trapping potential', which is nearly quadratic near the poles. When 
vortices conglomerate at the poles, they may coalesce to form giant or ring-like vortices, and 
these dynamics and phases are the topic of the present paper. 

Three contributions to the energy should be kept in mind to interpret the phase diagrams 
obtained in our calculations. First, to create a vortex, the kinetic energy of the associated 
supercurrent (on the 2D spherical surface) should be taken into account. This contribution 
increases when two vortices with parallel vorticity are placed near each other, so it acts as a 
repulsion between the vortices. Thus, it tends to favor splitting of the giant vortices. Second, 
to create a vortex, the order parameter needs to be suppressed over a region typically of 
the size of the coherence length. The energy cost associated with this turns out to favor a 
multiply quantized (giant) vortex over the corresponding $o-iiiultivortex state. The energy 
cost is relatively larger for a smaller sphere, since proportionally a larger fraction of the total 
order parameter needs to be suppressed. The balance between these two energy contributions 
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FIG. 9: Phase boundaries for spherical superconducting layers with different radius TZ and thickness 
W as a function of the temperature T and the applied magnetic field H. The boundaries between 
the thermodynamically stable normal state, Meissner state, single $o-vortex state, giant vortex 
state, and $o-™ultivortcx states arc shown by solid lines. The dashed line approximately indicates 
the boundary between the regions, where giant vortex states or $o-™ultivortex states are the 
thermodynamically stable state. 
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can be used to qualitatively understand the phase diagrams that we calculate for thin shells. 
Indeed, for magnetic fields corresponding to multiple quanta of vorticity, the smaller spheres 
will favor giant vortices, whereas the larger spheres favor the <l>o-niultivortex state. Note 
that this contribution to the energy strongly disfavors ring-like vortex states. 

The third contribution to the energy is related to the gradients in the magnetic field. 
When the shell is much thinner than the penetration depth, the currents on the shell will 
not substantially perturb the applied field, and this contribution plays no role. However, 
for thicker shells, this contribution does become important - as can be seen from Fig. [7], 
the magnetic field is substantially perturbed. When a $o-Kiultivortex lattice is present, the 
magnetic field flux is concentrated near each vortex core, and shielded in between, leading 
to a larger magnetic contribution to the energy than for a giant vortex. Thus, for a thick 
shell, this contribution will favor the giant vortex state. This agrees with our phase diagram 
showing that the region, where the giant vortex is stable, is growing for thicker shells. 

The temperature dependence of the phase diagrams was studied straightforwardly by 
taking temperature into account through the Ginzburg-Landau parameters. When multi- 
ple quanta of vorticity are present, we find that the giant vortex phase forms the preferred 
high-temperature phase. This offers the prospect of probing a temperature-driven transition 
between a giant vortex and a $o-niultivortex state, alongside with a magnetic-field driven 
transition. Moreover, the vortex dynamics are shown to be not sensitive to moderate im- 
perfections in the shell; the energy contributions discussed here can overcome the pinning 
potential due to for example thickness inhomogeneities - such pinning potentials have in 
past experimental work hampered the detection of the giant vortex state. This robustness, 
together with the tunability of the phase diagram through a limited set of controllable pa- 
rameters, makes superconducting nanoshells uniquely suited for the study of novel vortex 
states. 
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